Effects of axial torsion on sp carbon atomic nanowires 

Luca Ravagnan^'^ Nicola Manini^'^, Eugenio Cinquanta^'^*, 
Giovanni Onida^'^jH Davide Sangalli^'^, Carlo Motta^'^, Michele 
Devetta^'^, Andrea Bordoni^'^, Paolo Piseri^'^, and Paolo Milani^'^ 
^Physics Department and INFM - University of Milan, 

Via Celoria 16, 20133 Milano, Italy 
^CIMAINA - Via Celoria 16, 20133 Milano, Italy and 
^European Theoretical Spectroscopy Facility, 
Via Celoria 16, 20133 Milano, Italy 
(Dated: April 22, 2009) 

Abstract 

Ab-initio calculations within Density Functional Theory combined with experimental Raman 
spectra on cluster-beam deposited pure carbon films provide a consistent picture of sp-carbon 
chains stabilized by sp^ or sp^ terminations, the latter being sensitive to torsional strain. This 
unexplored effect promises many exciting applications since it allows one to modify the conductive 
states near the Fermi level and to switch on and off the on-chain 7r-electron magnetism. 
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Post-silicon electronics has seen the recent opening of entirely new perspectives along the 
way of carbon-based devices. By proper nanoscale design, entirely carbon-made transistors 
have been realized Future applications have been devised, including bio-nanotechnology 
ones such as devices for fast DNA reading 21]. Even considering only well-demonstrated 
applications, the potential of carbon-based electronics is undoubtedly enormous, as testifiec 
by the realization of non- volatile memories based on two-terminals atomic-scale switches [Si 

n 

and bistable graphitic memories [4| . Specifically, these structural memory effects have been 
explained by the formation of carbon chains made by a few aligned sp-hybridized atoms 
bridging a nanometric gap 

In this context, the production of pure carbon nanostructured films with coexisting sp 
and sp^ hybridization s], 0] opens the exciting possibility to tailor complex carbon-based 
nanostructures with linear chains made of sp-hybridized C atoms connecting graphitic nano- 
objects. However, despite several theoretical studies devoted to sp carbon nanowires 



10[ | classified either as cumulenes (virtually conducting, characterized by double C-C bonds) 



or polyynes (large-gap insulators with alternating single and triple bonds), the implications 
associated with the nanoscale geometrical manipulation of hybrid sp-|-sp^ carbon systems 
are still largely unexplored. 

In this Letter we show that sp nanowires can be stabilized effectively by termination 
on graphitic nanofragments, and that in the resulting structures the ~ 1 nm-long linear 
atomic chains can be torsionally stiff, due to the broken axial symmetry with staggered 
TT-bonds. This stiffness is rich of consequences. We explore here how the structural, vibra- 
tional and electronic properties of such chains are affected by the nature and geometry of the 
termination. In particular, we show that sp^ bonding to graphitic fragments and graphene 
nanoribbons (NRs) produces remarkably stable structures, with cumulene-type chains dis- 
playing a non-negligible bond length alternation (BLA), so that the traditional categories 
of polyynes (alternating single-triple bonds, yielding a large BLA) and cumulenes (double 
bonds, negligible BLA) appear too simplistic for the description of these systems. Torsional 
deformations are found to affect the BLA, electronic gap, stretching vibrational frequencies 
and spin magnetization of the chains. 

We study these effects in realistic nanostructures, including carbon chains bound to 
graphitic fragments. We perform all calculations within Density Functional Theory in 
the Local Spin Density Approximation, using a plane-wave basis as implemented in the 
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FIG. 1: (Color online) A few representative structures involving an 8-atoms sp-bonded carbon 
chain terminated on sp^ carbon fragments [(a-e): NRs; (f): C2o]- Either edge termination [(a-d), 
sp^-like] or termination on an internal atom of the fragment [(e-f), sp^-like] is possible. Binding 
energies (with respect to the uncapped straight chain plus fully relaxed sp^ fragments) and BLA 
are reported. 

ESPRESSO [llj suite ^4\. 

Figure [1] displays a few of the studied systems involving either sp^ or sp'^ bonding of a sp 
nanowire with an sp^-type fragment. The chosen end-capping nanostructures include planar 
graphitic fragments and closed-cages clusters (here, C20, the most curved fuUerene). These 
structures are intended to represent typical interfaces present in the nanostructured films 
produced by cluster-beam deposition |12|. 

The nature of the terminal bonding turns out to be crucial in determining the structure 
and electronic properties of the wire. An sp^-kind termination produces remarkably stable 
cumulene-type structures (between 7.9 and 13 eV for the formation of the two new bonds), 



characterized by a BLA between 5 and 9 pm [25|]. The computed binding energy should 
be compared with the energy-per-bond of 2.1 eV that we obtain for a lateral attachment 
of the same chain to the ribbon edge, and with the formation energy of graphene edges 
moreover, it is much larger than the reported binding energies of carbon chains 



inside nanotubes [15|. Figure [T] also shows that the mere value of the BLA does not allow 
one to distinguish between carbon chains which would be traditionally classified as cumulenic 
(a-d) or polyynic (e-f) according to their terminations. 

For the sake of comparison we also consider standard cumulenes and polyynes, in the form 
of isolated carbon chains stabilized by hydrogen terminations. Polyynes, which can be seen 
conceptually as acetylene molecules with longer alternating triple/single-bond carbon chains 
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(C„H2), have been synthesized up to a considerable length (n = 20) in liquid and solid 

matrices, and also with different stabilizing end-groups. Their electronic and vibrational 
properties as isolated species have been characterized extensively, mainly by electronic and 



Raman spectroscopy 



17 



18|. On the other hand, cumulenes C.„+2H4 can be seen as sp 



chains terminated by CH2 groups, yielding all double C=C bonds, and can be ideally thought 
as ethylene molecules lengthened by extra carbon atoms. Cumulenes are more elusive and 
less well characterized than polyynes, due to their fragility. Recently, short cumulenic chains 
have been synthesized in their basic forms, butatriene and hexapentaene [lol. In fact, 
cumulenic chains are often produced in conjunction with more complex terminations than 
simple CH2 units, including CPh2, i.e. 1,1 diphenyl ethyl (DPE) groups ^20|], which we also 
simulate. 

As one could infer from elementary valence bond or tight-binding considerations, depend- 
ing on the number n of carbons being even or odd, sp^-terminated cumulenes assume a D2/1 
(planar) or a Y)2d (staggered) ground-state geometry respectively, due to the alternating ori- 
entation of the TT bonds along the chain 19|]. Similarly, chains bonded to sp^ structures are 
affected by the relative orientation of their terminations. Indeed, a memory of the orienta- 
tion of the bonds of the terminating sp^ carbon propagates along the sp-hybridized chain, so 
that even-n chains tend to relax to a configuration where the termination sp^ planes coincide, 
while odd-n chains tend to keep their terminations at a twist angle 6 = 90°. As a conse- 
quence, despite their purely one-dimensional nature, sp^-terminated carbon chains display 
a non-vanishing torsional stiffness, no matter if they are straight or bent as in Fig. [TJ:;,d. In 
contrast, ideally polyynic chains (i.e. those terminated at a sp^ site, with a pure single-triple 
bond alternation) are almost completely free to rotate around their axis, but suffer from an 
obvious frustration when the number of atoms is odd since the long-short bond alternation 
must swap at their middle |2l|. Importantly, in nanostructured cluster-assembled carbon 
characterized by a complex three-dimensional arrangement of graphitic fragments and sp 



chains j^, Q], a large number of the chains binding to sp^ structures are not free to relax 
their terminations to the preferred angular geometry, and must hence be expected to be, in 
general, strained torsionally. 

Since the simulation of carbon-only structures such as graphene NR bridged by chains 
allow us to investigate few relative angular arrangements only, we extend our study also 
to chains with simpler saturating ligands, namely CII2 and DPE. The latter turns out to 
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FIG. 2: (Color online) Total torsional energy (a), Kohn-Sham electronic gap (b), and bond-length 
alternation (c) as a function of the twist angle 9 for representative even-numbered sp-carbon chains 
with different terminations. Open/filled symbols refer to the low/high-spin electronic configura- 
tions. 

reproduce better the behavior of a large (potentially semi-infinite) graphitic fragment, which, 
at variance with CH2, shares with the chain only a fraction of its unsaturated p^ electron, 
which is partly delocalized over an extended aromatic sp^ structure. 

Figure [2] summarizes the influence of different end groups on the BLA and, for sp^ termi- 
nation, the torsional strain energy and the Kohn-Sham electronic gap of even-n chains as a 
function of 9. Interestingly, largely strained chains undergo a magnetic instability, turning 
spin polarized. The reason is the quasi-degeneracy of two yr-bonding/antibonding electronic 
levels near the Fermi energy illustrated by the closing of the gap. Fig. [2)3. Remarkably, 
in even-numbered chains of all considered lengths, no matter if CH2- or DPE-terminated, 
the switching to a spin-polarized configuration takes place at the same twist angle Og — 79°, 
highlighted in Fig. [2] [26!]. This Og invariance implies that the energy gap A and the exchange 
splitting of the electrons near the Fermi level scale in the same way. 

Calculations show that the BLA of the sp^-terminated chains varies substantially with 
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the nature of the termination itself 



10|. The length of the extremal bond (i.e. the one 



connecting the last sp carbon with sp^/sp^-hybridized ligand), which correlates with the 
BLA, is minimal in the case of a simple CH2 termination, but increases substantially in DPE- 
terminated chains, assumes even larger values in NR-terminated wires, and is maximum for 
polyynic-type terminations, see Fig. [T] and Fig. The torsional barrier is consistently 
smaller for NR terminations, as indicated by the cross at 6' = 90° in Fig. [2^. 

In the light of the above results, chain-termination details are expected to influence the 
vibrational properties as well. The latter offer an invaluable opportunity to check if the the 
considered structures are representative of those present in nanostructured cluster-assembled 
films, for which Raman spectra are the main experimental evidence of the presence of linear 
carbon chains. In fact, in previous works some of us showed that the Raman fingerprint of 
carbyne chains in sp-sp^ carbon is characterized by a broad feature, where 2 components CI 
and C2 peaked at 1980 and 2100 cm~^ respectively can be recognized . Traditionally 
these features were attributed generically to cumulenes (CI) and polyynes (C2). We hence 
calculate the phonon frequencies and eigenvectors of the structures exemplified in Figs. [1^ 
and [Tf, plus CH2- and DPI5-terminated carbynes of several lengths, using standard Density- 
Functional Perturbation Theory [llij^^]. As a benchmark, our theoretical C-C streching 



modes of polyynes C„H2 (n = 8 — 12) match the experimental frequencies [18|| to within 
40 cm"^ 

Beside several bending and long-wavelength stretching modes, whose low frequency falls 
in the same range as the vibrations of graphitic and diamond-like carbon material, short 
linear carbon chains display a few characteristic Raman-active stretching modes in the range 
1800 -T- 2300 cm~^. One mode, sometimes named the a-mode in the literature 



shows 



a displacement pattern localized near the chain center, and usually bears the strongest 
Raman intensity Since the displacements at the chain ends are less than 10% of those 
of the central atoms, the frequency of the a-mode is almost unaffected by the mass of the 
termination (e.g. calculations for C6H4 and C8H4 with 1000-times increased hydrogen mass 
give frequency shifts of less than 1 cm~^). 

The stretching frequencies of sp chains turn out to be influenced by: (i) the type of 
termination (sp^ vs. sp^); (ii) the chain length, with even/odd alternation effects; (iii) for sp^ 
termination, the relative orientation of the termination themselves, with effects of torsional 
strain. The calculated frequencies of the high-energy Raman-active modes display a clear 
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FIG. 3: (Color online) (a) The computed frequency of the Raman a-mode (horizontal scale) 
for carbon chains of different length n (vertical scale) and with different terminations (b) The 
softening of this mode for CH2-terminated chains as a function of the twist angle (vertical scale) 
(c) The experimental Raman spectra of pristine cluster-assembled sp-sp^ film (grey line) and of 
the same material after 2 days exposure to He, 100 Torr (black line). The underlying Gaussians 
report the empirical analysis of both spectra, resulting in 5 components at frequencies separated 
by approximately 80 cm~^. The individual components display different decays, beside becoming 
narrower and undergoing a ~ 10 cm^^ blueshift. 

distinction between sp^ and sp'^-terminated chains, as shown in Fig. [31 Only even-numbered 
chains are reported, since odd chains lack a modes, and have in general much smaller Raman 



cross sections 



28|. 



Figure [3b displays the in-situ Raman spectrum of an sp-sp^ nanostructured-carbon (ns-C) 
film '5| in the carbyne region, measured using the 488 nm line of an Ar"*" laser and acquired 
with very high signal-to-noise ratio. The spectrum of the as-deposited material is compared 
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to that obtained after exposure of the film to He in order to promote sp chain decay 22 1. 
Clearly, a description in terms of two peaks only cannot account for the complex structure 
and decay observed. In particular, the Cp component at the highest frequency (peaked at 
2260 cm~^) can be attributed uniquely to short polyynic chains, as it is higher than any 
cumulenic a mode (see Fig.[2t'), while the other components can be related both to polyynes 
and cumulenes of different chain length. As illustrated in Fig. after He exposal, individual 
components have different evolutions during the C band decay, and in particular the peaks 
at lower energy, corresponding to longer chains, decay faster than the higher-energy ones 
(i.e. shorter chains). Furthermore, the Cp peak does not shift during the decay nor change 
its width, while all lower peaks are blue-shifted by ~ 10 cm^^ and narrowed by ~ 7 cm^^. 
Indeed calculations, summarized in Fig. ^jp, show that the high-frequency stretching modes 
of torsionally-strained CH2-terminated chains are affected quite strongly by the twist angle, 
with a redshift up to ~ 100 cm~^. However, since chains with smaller torsional barrier 
(such as those bound to DPE and nanoribbons) show smaller redshifts, this effect evaluated 
for CH2-terminated chains should be considered as an upper limit for realistic pure-carbon 
nanostructures. 

The observed blue shift of the peaks accompanying the decay can then be explained 
if each peak is related to a particular family of cumulenes, having all the same length but 
different strain: the more strained chains, having softer Raman modes, decay faster than the 
others, resulting in a net blue shift and narrowing of the peak. A faster decay of torsionally 
strained vibrationally red-shifted cumulene-type is indeed to be expected due to their higher 
total energy (Fig. [2^). On the contrary, no torsional strain applies to polyynes, and this is 
why the Cp peak does not shift. 

In summary, we performed ab-initio total-energy and phonon calculations on a selected 
range of model structures sampling significantly the infinite variety of three-dimensional 
arrangements of linear carbon chains bridging graphitic fragments in different hybridization 
states. Theoretical results suggest that sp-carbon chains are stabilized in particular by 
bonding to the edges of graphitic nanofragments, and allow us to interpret the nontrivial 
features and decay of experimental Raman spectra of cluster-beam deposited pure carbon 
films. Moreover, the data for sp^-terminated chains point towards a rich phenomenology 
driven by even/odd alternation effects and by the effects of torsional strain. The latter 
modifies the electronic states near the Fermi level, suggesting the possibility to control the 
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nanowire conductance 



211 ]. optical properties, and spin magnetization, purely by twisting its 



sp^ termination, e.g. by coupling terminating graphene sheets with micromachined torsional 
actuators 23|. Linear carbon chains bridging graphene nanogaps, recently proposed as an 



explanation of the conductance switching in two-terminals graphene devices |3|, could 
hence acquire an important role in future carbon-based electronics. 

We acknowledge support by the Italian MIUR through PRIN-2006025747 and by the 
European Union through the ETSF-I3 project (Grant Agreement No. 211956 / ETSF User 
Project No. 62). 
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Ill ] with a wavefunc- 



[24] We use the C.pz-rrkjus.UPF ultrasoft pseudopotential from Ref. 

tion/charge cutoff of 30/240 Ry, and relax all atomic positions until the largest residual 
force is < 2 X 10"^ Ry/oo (8 pN). In the case of periodic graphene NRs, we adopt supercells 
with three hexagonal units along the periodic direction and at least 7 A of vacuum separating 
periodic images in the two other directions, optimizing the lattice constant until the stress 
tensor drops below 2 x 10^^ Ry/af]. We sample the Brillouin zone with at least 13 A;-points 
in each periodic direction and only /c = in non-periodic directions. Numerical details are 
similar to those validated and used, e.g., in Ref. [14]. 

[25] The BLA measures the degree of dimerization and, excluding the terminal bonds, can be 
defined as i Z]"=i(d2i-i + dn-{2j-i)) /^e - I]"=i(f^2i + dn~2j)/no , with di = |fj-r;+i|, Ue = 
(n + 2)/4, and rio = n/4 (integer part). 

[26] Odd-numbered chains show a reversed behavior, with high-spin states near the energetically 
unfavorable planar geometry ~ 0°. 

[27] All phonon calculations start by fully relaxing all degrees of freedom except, in the case of 
torsionally strained C„(CH2)2 systems, the angular coordinates of the H atoms around the 
molecular axis. 

[28] For the simplest structures (CH2-terminated chains) we also compute Raman intensities. 
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